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FOREWORD 


This report presents results of theoretical research on Atmospheric 
Inferences from Infrared Horizon Measurements performed under National 
Aeronautics and Space Administration Contract NAS1-7572 for Langley 
Research Center. 

This study is concerned with the problem of extracting information on 
atmospheric structure from infrared radiance measurements near the earth f s 
horizon. The study includes an assessment of the 15(i CO 2 band and the H 2 O 
rotational band as potential spectral regions whose radiance measurements 
can be used to infer temperature and moisture profiles in the stratosphere. 
Investigations including inversion techniques, absorption band structures, 
error analyses and problems relating to radiometer design constitute por- 
tions of the study. 

Gratitude is extended to the scientists at NASA/Langley Research Cen- 
ter for their interest, technical assistance and guidance during the course 
of the study. These scientists include Mr. R. Davis who is technical mon- 
itor of the program and official government representative concerning any 
questions and availability of the results, and also include Mrs. R. Whitman 
and Messrs. T. McKee, J. Dodgen and H. Curfman. 
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ABSTRACT 


This study concerns the feasibility of inferring temperature and mois- 
ture profiles in the stratosphere from horizon radiances in the 15 micron 
CO 2 band and the H 2 O rotational band. Investigations are performed of in- 
version techniques, absorption band structures, error analyses and prob- 
lems related to instrument design. A relatively simple inference method 
termed the "instant inversion" technique is developed and applied success- 
fully to both theoretical radiance profiles and observed horizon radiances 
from 'Project Scanner. A detailed analysis of these inversion results leads 
to a recommended measurement and data reduction technique to determine both 
temperature and moisture profiles in the stratosphere and lower mesosphere. 
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INFERENCE OF STRATOSPHERIC TEMPERATURE AND MOISTURE PROFILES 
FROM OBSERVATIONS OF THE INFRARED HORIZON 


By Frederick B. House and George Ohring 
GCA Corporation, GCA Technology Division 
Bedford, Massachusetts 


SUMMARY 


Infrared radiance emitted by the atmosphere above the earth's hard 
horizon varies with tangent height depending on the concentration and 
spectral absorbing properties of the atmospheric gases and on their pres- 
sure and temperature distribution with altitude. Therefore, measured hori- 
zon radiances contain information about these atmospheric parameters. This 
study concerns the feasibility of inferring temperature and moisture pro- 
files from horizon radiances in the 1 5(a CO 2 band and the H 2 O rotational 
band . 


Investigations are performed of inversion techniques, absorption band 
structures, error analyses and problems related to instrument design. A 
relatively simple method termed the "instant inversion" technique is de- 
veloped and tested extensively. The successful application of this tech- 
nique to both theoretical profiles and observed horizon radiances from 
Project Scanner demonstrates the feasibility of inferring atmospheric 
parameters. A detailed analysis of these inversion results leads to a 
recommended measurement and data reduction technique to determine both 
temperature and moisture profiles in the stratosphere and lower mesosphere. 
It is concluded that useful temperatures can be obtained in the altitude 
range from about 10 to 65 km and water vapor mixing ratios can be deter- 
mined within the 10 to 50 km altitude range. 
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INTRODUCTION 


The inadequacy of present stratospheric measurements of temperature 
and moisture are apparent when one considers available instrumentation, 
geographical coverage and observational platforms for in-situ measurements. 
Temperature profiles are available in the lower stratosphere up to alti- 
tudes of 25 to 30 km from standard meteorological balloon soundings. Vari- 
ous rocket-borne instruments have been developed for sounding the atmo- 
sphere between 30 to 55 km, with occasional measurements to higher alti- 
tudes. Moisture profiles are made infrequently in the stratosphere, and 
measurements require sophisticated instruments such as frost-point hygro- 
meters. When one considers the frequency and geographical coverage of 
these measurements on a world-wide basis, the lack of routine observations 
over oceanic regions and in the Southern Hemisphere is quite pronounced, 
even in the troposphere. 

The need for a better specification of the general state of the at- 
mosphere at any time is indicated by numerical experiments using general 
circulation models (ref. 1 and 2). Prediction models currently under de- 
velopment require temperature, moisture and wind information up to alti- 
tudes of 40 km (ref. 3). The requirement for data in the upper atmosphere 
will increase as prediction models become more sophisticated. 

Possibilities of remotely probing the atmosphere from satellite plat- 
forms are attractive in view of the requirement of global coverage of mete- 
orological data. The second generation of meteorological satellites in 
the NIMBUS series will incorporate sophisticated radiometers for remote 
probing of the atmosphere in the infrared region of the electromagnetic 
spectrum. Outgoing radiance signals from atmospheric absorption bands de- 
pend on the concentration and spectral absorbing properties of the gases 
and on the pressure and temperature distribution with altitude. There- 
fore, radiance measurements contain information about these atmospheric 
parameters. The inference of parameters from radiance measurements is 
commonly termed "the inversion problem." 

Active absorption bands in the infrared region of the spectrum in- 
clude carbon dioxide (CO 2 ) , water vapor (H 2 O) and ozone (O 3 ) . Assuming 
CO 2 is uniformly mixed in the atmosphere, it is possible to infer verti- 
cal temperature profiles from radiance measurements in the 1 5p CO 2 band. 
With knowledge of the temperature profile, it is then possible to deter- 
mine vertical concentration profiles of water vapor and ozone based on 
radiance measurements in the 6.3 or 20 to 40p bands of H 2 O and in the 9 . 6 \± 
O 3 band, respectively. 
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Measurements by current nadir-viewing radiometers will allow accurate 
temperature inferences to altitudes from the surface to about 30 km (ref, 
4). Moisture inferences are possible only in the troposphere (ref. 5), 
and determination of total O 3 concentration and an estimate of its alti- 
tude distribution (ref. 6 ) appear feasible from nadir experiments. 

The possibility of inferring atmospheric parameters from horizon radi 
ance measurements has not received the attention given to nadir-viewing ex 
periments. The reason for this probably lies in the fact that temperature 
and moisture data in the stratosphere are not as useful as data in the 
troposphere, and that inferences apply to the region of the atmosphere 
above cloud-top altitudes. It should be pointed out, however, that tropo- 
spheric inferences of temperature suffer from the same cloud problem (ref. 
7). Inferences of atmospheric parameters from horizon radiances have sev- 
eral important advantages: excellent vertical resolution, enhanced signal 

due to the long optical paths associated with limb viewing, and no back- 
ground interference since the radiances are viewed against the cold of 
space. Data derived from horizon radiances certainly would be a compli- 
mentary adjunct to current nadir-viewing experiments. 

The current study constitutes a detailed investigation of the temper- 
ature and moisture inference problem from observations of the infrared 
horizon. Work along similar lines has been performed concurrently at the 
Langley Research Center (ref. 8 and 9), and the possibility of the concept 
was apparent in an investigation of horizon definition (ref. 10 ) for space 
craft attitude control. The study includes an assessment of the 1 5p CO 2 
band and the H 2 O rotational band as potential spectral regions whose radi- 
ance measurements can be used to infer temperature and moisture profiles 
in the stratosphere. Investigations of inversion techniques, absorption 
band structures, error analyses and problems relating to radiometer de- 
sign constitute parts of the study. 

Benefits to be realized from routine measurements of temperature, 
moisture and ozone profiles in the stratosphere are many. These data 
would aid research on the nature of stratospheric warmings, troposphere- 
stratosphere interactions, the 26-month cycle in the equatorial strato- 
spheric winds, diurnal variation at the stratopause (ref. 11 ), strato- 
spheric general circulation, global distribution of moisture and ozone, 
radiative balance in the stratosphere and lower mesosphere, circulation 
of trace constituents and the photochemistry of ozone. Practical appli- 
cations of these data include operational problems associated with high- 
flying aircraft like the SST, effect on development of numerical predic- 
tion models and problems related to National Defense including atmospheric 
effects on the operation of offensive and defensive missile systems. 
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SYMBOLS 


constant in an equation for a straight line (equation 15) 

A^ slope coefficient in an equation for a straight line (equation 15) 

2 -1 

B Planck radiance, watts/m -sr-cm 

— 2 
B effective Planck radiance for a given tangent height, watts/m -sr 

•“8 4 2 

constant in Planck formula, 1.1909 x 10 watts-cm /m -sr 
constant in Planck formula, 1.4389 cm deg 
h tangent height, km 

i counting index 

2 -1 

J source function, watts/m -sr-cm 

L generalized absorption coefficient 

2 

N horizon radiance, watts/m -sr 

n counting index 

P pressure, mbar 

re random error 

rms root-mean-square error 

se systematic error 

T temperature, °K 

— o 

T effective emitting temperature, K 

* 

U reduced optical path length of a gas 

W mixing ratio of a gas, g/kg 

WF weighting function, k ^ 

x independent variable along a ray path 

z altitude variable, km 

z' altitude of the effective emitting level, km 
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A increment of a variable 

e emittance 

v wavenumber, cm ^ 

V c central wavenumber of a bandpass filter, 

v ^ lower wavenumber of a spectral interval, 

upper wavenumber of a spectral interval, 



T transmittance 


Subscripts 

a anterior portion of ray path 

C carbon dioxide 

H water vapor 

p posterior portion of ray path 

v at a wavenumber 

Av specified spectral interval 


-1 

cm 

-1 

cm 

-1 

cm 
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INSTANT INVERSION TECHNIQUE 


Background 

A distinction should be made between the concepts of conventional in- 
situ measurement and remote sensing of atmospheric parameters. In-situ 
measurements of atmospheric parameters such as temperature, pressure and 
moisture are usually made by separate instruments which allow a point -by- 
point determination of each parameter as a function of altitude. The in- 
struments are carried through the atmosphere by some platform such as a 
balloon, parachute or rocket. The concept of remote atmospheric sensing 
is quite different in that the vertical profile of radiant energy, emitted 
by the atmosphere, is the only measured quantity. A single observation of 
radiance in a given absorption band of an atmospheric gas is related to a 
conglomeration of many parameters at different altitudes of the atmosphere. 
The magnitude of a radiance measurement is related principally to temper- 
ature, pressure, concentration of absorbing gas and the optical properties 
of the absorption band. There is indeed a difference between the concepts 
of in-situ measurement and remote radiometric sensing of atmospheric 
parameters . 

As stated earlier, the objective of this study is the extraction of 
two interesting parameters, temperature and moisture, from radiance obser- 
vations of the infrared horizon. In order to determine either of these 
parameters, it is necessary to have knowledge of all other parameters af- 
fecting the observed radiance. For example, to determine temperature it 
is necessary to have knowledge about pressure, the amount of gas and its 
optical properties. If the determination is of moisture or the concentra- 
tion of water vapor, it is necessary to know the distribution of pressure 
and temperature in the atmosphere as well as water vapor's optical proper- 
ties. The whole problem is further complicated by the fact that a single 
radiance observation is an integral of the emitted energy from many levels 
in the atmosphere. 

The spectral regions of interest in this study include the CO 2 band 
from 615 to 715 cm“l and the H 2 O band from 315 to 475 cm”l. Examples of 
radiance profiles in these bands are presented in figure 1. 

The determination of atmospheric parameters from radiance measure- 
ments is commonly termed the inversion problem. The discussion of this 
section presents the theoretical background and justification of a simple 
method termed the "instant inversion technique." This technique is de- 
veloped for both temperature and moisture inferences, and it takes advan- 
tage of the limb-viewing geometry associated with horizon radiance measure- 
ments in the infrared region of the electromagnetic spectrum. 
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HORIZON RADIANCE N (hi. Wolfs/ m sr 



Figure 1. Measured horizon radiance profiles in the CO 2 band (615 to 
715 cm“l) and the H 2 O band (315 to 475 cm“l) from Project 
Scanner observations (ref. 8). 
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Temperature Inversion Theory 

The horizon radiance that would be sensed by a satellite-borne 
radiometer in figure 2 is given by the exact form of the radiative trans- 
fer equation 


2 V h >-’ 

(h) = '/ f V V T(X)] d %( h > x ) dv (1) 


where it is assumed that the atmosphere is the only source of emission 
and that there are no obstructions such as clouds or the earth's surface. 
Equation (1) indicates that the horizon radiance is the sum of incremental 
energy sources at temperature T integrated long the line -of -sight, each 
having a change in spectral transmittance of dT y . The measured radiance 
arising from a given absorption band is integrated over the entire band 
between limits and • The inversion problem consists in the extrac- 
tion of temperature information from the integral equation (1). 

As long as the assumption of local thermodynamic equilibrium is 
valid, the source function J v [T(z)] is the Planck function 

C 1 v 3 

J [T(z) ] - B [T (z) ] = (2) 

v s v L N / CL v \ v 

exp -4- - 1 


Thus, after integration over the band, equation (1) may be expressed by 

%<*■> ■ J B Av [T(z)1 ST dx < 3 > 

+00 

It is assumed in the instant inversion technique that the Planck radi- 
ance B^ v may be expressed by an average Planck radiance B^ v [T(z')], associ- 
ated with a single temperature at a specified altitude z'. With this 
approximation, equation (3) is easily integrated 


N Av (h > * • T Av> B Av tI(2 ' )5 ’ e Av (h) B Avt T(2 '» 

Solving for the Planck radiance in equation (4) gives 


N (h) 


( 4 ) 


( 5 ) 
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which is interpreted in terms of temperature using a form of equation (2) 
after including the spectral band width. 

In this study the coordinate x does not extend from the dimensions 
+oo to but is bounded by a spherical shell at an altitude of 90 km 
above the earth's surface. The total path length emittance e^ v (h) in 
equation (5) may be expressed as the sum of two integrals representing 
the contributions by the portions of the atmosphere anterior and posterior 
to the tangent height point of the ray path. 

6 Av (h > - • f(-sr) dz - f (-sr) dz (6> 

90 a h P 


The weighting function or change of transmittance along the ray path 
with respect to the altitude coordinate is of particular interest in de- 
picting the physics of the problem. Assuming spherical symmetry, the 
summed weighting function WF(h,z) at a given altitude is the sum of an- 
terior and posterior contributions 

-fe) < 7 > 

a p 

Figure 3 illustrates the variation in weighting functions with altitude 
for different tangent heights. The curves give a good approximation of 
the contribution by different levels of the atmosphere to the radiance 
at a given tangent height. As can be seen by the peaked nature of the 
functions, the radiation for tangent heights above 25 km originates from 
a rather narrow altitude increment above the reference level. This 
optical property of the atmosphere is unique to limb-viewing and justi- 
fies the basic assumption of the instant inversion technique (equation 5), 
that the Planck radiance for a representative temperature can be inferred 
from the ratio of the observed radiance and the total emittance. 


Moisture Inversion Theory 


The approach to the inference of moisture profiles follows a similar 
line of thought as that for temperature inferences. However, the total 
emittance of water vapor is the computed quantity rather than the Planck 
radiance, and the mixing ratio is one parameter in an exponential 
expression for the total emittance. The emittance is computed by 


N H (h) 

e H (h) = B h [T(z')] 


( 8 ) 
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ALTITUDE 




where the subscript (H) refers to the rotational water vapor band (315 to 
475 cm~l) and the Planck radiance is determined from the temperature 
derived previously from the CO 2 inversion. In a nonscattering atmosphere, 
the transmit tances of l^O and CO 2 are related to the emittances by 


and 


T H (h) = 1 " e H (h) 

(9) 

T c (h) = 1 " £ c (h) 

(10) 


where the subscript (C) refers to the carbon dioxide band (615 to 715 cm”^). 

The problem at hand is to establish the relationships between the two 
transmittances in (9) and (10) based on transmission theory for the respec- 
tive absorption bands. A convenient transmission model for this discussion 
is that developed by Elsasser (ref. 11) where the transmittances are ex- 
pressed by the sum of common logarithms, Log U* + Log L. In this case U* 
is the reduced optical path length of the gas and L is the generalized 
absorption coefficient. The path lengths for each gas are defined as 



where the capital sigma is a shortened notation for the summation term 
across n shells along the line of sight, and P Q and T 0 are standard 
pressure and temperature, respectively. It is noted that the summation 
term is the same for both gases when it is assumed that the mixing ratio 
of water vapor is constant above the tangent height altitude. Units of 
optical path length are cm for CO^ and grn/cm^ (or precipitable cm) for 
H 2 O. Since the total emittance of CO 2 is known, it is possible to 

evaluate the term ^ in equation (11) 


U 

W~ 


(13) 
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when it is assumed that CO 2 is uniformly mixed in the stratosphere. Sub- 
stitution of equation (13) in (12) and solving for the mixing ratio gives 


U W 

H C 

1.292 x 10“ 3 U* 


(14) 


Both of the optical path lengths U H and U q are determined from their 
respective emittances with knowledge of the band-pass filters of the radi- 
ometers . 


The above technique for inferring moisture profiles will work effec- 
tively provided that the H 2 O and CO 2 band-pass filters on the radiometer 
are matched optically. This statement implies that the optical path 
lengths for both gases are approximately the same as a function of alti- 
tude. Putting this condition in another context, this statement implies 
that the summed weighting functions in figure 3 peak at the same alti- 
tudes and that the radiance curves in figure 1 should increase in magni- 
tude at about the same altitudes. 

The CO 2 band (613 to 715 cm~l) and the H 2 O band (315 to 475 cm - !) 
measured by Project Scanner radiometers are not optically suited for water 
vapor inferences using the solution described by equation (14). However, 
the moisture inference concept may be investigated by utilizing temper- 
ature and pressure profiles based on meteorological soundings of the at- 
mosphere at the time of observation. These data allow a direct computa- 
tion of the U*c/Wc term in equation (14) rather than deriving this term 
from the total emittance parameter of C 02 * 

Summary Remarks 

In the instant inversion technique the assumption is made that the 
radiative transfer equation may be approximated by the product of the 
total path length emittance and the average Planck radiance. There is 
physical justification for this assumption since the weighting functions 
for the CO 2 band indicate that the emitted radiant power originates from 
a rather narrow altitude region above the tangent height level. With 
this approximation it is possible to determine atmospheric temperature 
from the ratio of observed radiance to emittance. The technique may then 
be applied to radiance measurements in the water vapor rotational band to 
infer moisture profiles once the temperature structure is known providing 
that the two bands are matched optically. 
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TECHNIQUE DEVELOPMENT 


In this section a number of physical aspects of the inversion prob- 
lem are considered which relate to the transmittance properties of the 
atmosphere in the CO 2 and H 2 O rotational bands. The total emittance in 
the instant inversion technique depends on the concentration of gas, the 
pressure and temperature or density along the ray path and the optical 
characteristics of each band. The value of emittance at a given level 
is influenced by the atmospheric structure above the tangent height level. 
Of course the H 2 O emittance varies with the magnitude of water vapor mix- 
ing ratio. Consideration is also given in this section to the dispersion 
characteristics of the horizon radiance with altitude and the effects of 
the lack of thermodynamic equilibrium. 


Altitude Dispersion 

The technique of temperature inversion should work quite well for 
all tangent heights above 20 km because the weighting functions are def- 
initely peaked above this level. Consider the curves of figure 4 which 
illustrate the same information as is in the weighting functions (figure 
3) but in a different manner. The altitude dispersion or percent of total 
radiance and its relation to the altitude above the tangent height level 
is shown by the curves. For example, the 50-percent radiance level for 
a 20-km tangent height is about 10 km (at a z' height at 30 km). In sharp 
contrast, the 50-percent radiance level for a 60-km tangent height is about 
1.3 km. The curves of figure 4 also illustrate the fact that the shape of 
the weighting functions is constantly changing with altitude. In prin- 
ciple, the physics of the atmosphere are such that the higher the altitude 
the better the temperature inference using the instant inversion technique. 

In all the results compiled in this study, the 50-percent radiance 
level for each tangent height is assumed to be the effective emitting 
level (z') for purposes of the instant inversion technique [see equation 

( 5 )]. 


Emittance Calibration 

The total emittance in the CO 2 band depends basically on the temper- 
ature and pressure structure of the atmosphere. Knowledge of the emittance 
variation with tangent height for a given radiance profile is essential to 
the success of the inversion technique [see equation (5)]. An emittance 
calibration technique was developed as part of this study to account for 
the variation in emittance with air mass structure. The technique depends 
on an additional piece of information which is assumed to be characteristic 
of a particular atmosphere. Two parameters were investigated for this pur- 
pose, the height of the 10 mbar surface and the summed radiances between 


14 



ALTITUDE ABOVE TANGENT HEIGHT Az(h), km 
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0 20 40 60 80 100 

PERCENT OF TOTAL RADIANCE 


Figure 4. Altitude dispersion of the horizon radiance for different^ 
tangent heights — ICAO Standard Atmosphere 615 to 715 cm 
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tangent heights of 20 to 60 km in steps of 5 km. Linear correlations 
were developed between the variation in emittance with tangent height for 
twenty climatological profiles (ref. 10) and the two calibration param- 
eters. Calculations of the emittance for these profiles were performed 
by the LRC technical monitor, Mr. R. E. Davis. 

Figure 5 illustrates the correlation of emittance with the height of 
the 10-mbar surface for a tangent height of 30 km. The plotted points 
are computed emittances for each of the 20 climatological profiles. A 
strong coupling between the emittance and the height of the 10-mbar sur- 
face is quite evident in the graph. The correlation between emittances 
and the 10-mbar surface for tangent heights from 20 to 60 km shows a sim- 
ilar linear coupling. Based on these data it is possible to account for 
the variation of emittances using an empirical relationship based on a 
straight line curve fit to the data at each tangent height. Thus, the 
emittance at a tangent height may be expressed by the equation 

e(h) = A- HT- A . - + A (15) 

1 10-mbar surface o 

A least squares fit to the data in figure 5 leads to computed constants 
of Ai = 0.051 km"l and A 0 = -0.930 when the height of the 10-mbar surface 
is given in kilometers. 

The results of linear curve fits to emittance data for different tan- 
gent heights are presented in Tables I and II where the emittance calibra- 
tion parameters are the height of the 10-mbar surface and the area under 
the radiance-tangent curve, respectively. Included in the tabulated data 
are the coefficients for equation (15), the statistical standard error of 
the estimate and the standard error percentage. 


Lack of Thermodynamic Equilibrium 

Lack of thermodynamic equilibrium implies that the source function 
in the radiative transfer equation is different from the Planck function 
[see equation (2)]. This problem was investigated using the model des- 
cribed in reference 10. Radiance profile calculations were made for the 
Standard Atmosphere, assuming two conditions; one where the lack of equi- 
librium conditions existed and the second where the Planck formula was 
assumed to be valid. The change in the magnitudes of radiances at each 
tangent height indicates the lack of equilibrium. The results are pre- 
sented in figure 6 where the ordinate values indicate the decrease in 
computed radiances due to the lack of thermodynamic equilibrium. The ef- 
fect is quite small in the altitude range of interest for the instant in- 
version technique, no more than 1.1 percent for the radiance at a 60-km 
tangent height. Correction for this effect is accomplished by dividing 
the observed horizon radiances N^ v (h) in equation (5) by the ordinate 
values for each tangent height, represented by the curve in figure 6. 


16 



TOTAL EMITTANCE AT 30kmC(b =30) 


I 



Figure 5. Correlation of total emittances with heights of the 10 mbar 
surface for a tangent height of 30 km — 615 to 715 cm"^. 
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TABLE I. - RESULTS OF LINEAR CURVE FITS TO EMITTANCE 
DATA FOR DIFFERENT TANGENT HEIGHTS 

Parameter - 10 mbar Surface 


Tangent 
Height , km 

Curve Fit 
Slope, A^ 

Coefficients 
* Const, Aq 

Standard Error 
of Estimate 

Standard Error 
Percentage 

20 

1.65-2 

4.35-1 

2.46-3 

0.26 

25 

3.50-2 

-2.52-1 

2.43-3 

0.29 

30 

5.10-2 

-9.30-1 

3.83-3 

0.59 

35 

5.48-2 

-1.25 

5.36-3 

1.21 

40 

4.80-2 

-1.20 

5.89-3 

2.10 

45 

3.54-2 

-9.27-1 

4.93-3 

2.92 

50 

2.30-2 

-6.13-1 

3.43-3 

3.43 

55 

1.40-2 

-3.76-1 

2.05-3 

3.56 

60 

7.48-3 

-1.99-1 

1.15-3 

3.56 


*Units of km _ l 


TABLE II. - RESULTS OF LINEAR CURVE FITS TO EMITTANCE 
DATA FOR DIFFERENT TANGENT HEIGHTS 

Parameter - Summed Radiances 


Tangent 
Height , km 

Curve Fit Coefficients 
Slope, A^* Const, A Q 

Standard Error 
of Estimate 

Standard Error 
Percentage 

20 

2.53-3 

8.87-1 

4.26-3 

0.45 

25 

5.15-3 

7 . 07 -1 

5.85-3 

0.71 

30 

7.64-3 

4.60-1 

6.92-3 

1.07 

35 

8.34-3 

2.38-1 

5.74-3 

1.29 

40 

7.40-3 

9.68-2 

3.84-3 

1.37 

45 

5.51-3 

3.25-2 

2.50-3 

1.48 

50 

3.61-3 

1.06-2 

1.26-3 

1.26 

55 

2.19-3 

3.36-3 

7.50-4 

1.30 

60 

1.17-3 

3.41-3 

5.57-4 

1.72 


*Units of (watts /m^-sr)"! 
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LACK OF THERMODYNAMIC EQUILIBRIUM 


1.000 


0.998 


0.996 


0394 


0.992 


0.990 


0.988 




ANGENT HEIGHT h, km 


dynamic equilibrium as a function of tangent 
Standard Atmosphere. 
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Comparison of CO^ and H^O Band Emittances 

An analysis was made of the band emittance properties for both CO 2 
and H 2 O in the spectral regions associated with the Project Scanner obser- 
vations. The transmission model of Elsasser (ref. 12) is convenient for 
purposes of this discussion. In this model the transmittances are expressed 
by the sum of common logarithms, Log U* + Log L, where lT is the reduced 
optical path length and L is a generalized absorption coefficient which 
varies in magnitude within each band. The band emittances were computed 
for both CO 2 and H 2 O — the results are presented in figure 7 for an effec- 
tive temperature of 233°K. The band emittances are plotted as a function 

of Log^ ^ term which is a parameter common to both gases and defined in 

equations (11) and (12). The family of curves for H 2 O illustrate the change 
in band emittance for different concentrations of water vapor mixing ratios. 

In addition to the results in figure 7, computations were completed 
for a number of different effective temperatures corresponding to those 
found in the stratosphere. It is possible to determine the sensitivity of 
band emittance to temperature from these data. Results indicated that the 
CO 2 band emittance was insensitive to temperature variations; however, the 
H 2 O band emittance exhibited a systematic shift of significant magnitude. 

The emittance increases with temperature for a given Log.^ ^ term in 

figure 7. An empirical formula was developed for the 1^0 band which ac- 
counts for this temperature effect. 

Log 10 ^ = Log 10 ^ + 0 . 0063 (T - 233) (16) 

T T=233°K 


The effect of the formula is to shift the family of curves in figure 7 
either to the left or right depending on whether the temperature is warmer 
or colder than 233°K, respectively. 

It should be noted that the spectral intervals for CCU and H^O are not 
matched optically for purposes of inferring mixing ratios from horizon 
radiances. A comparison of curves in figure 7 will illustrate this point. 

The rise in the CC^ band emittance occurs in the Log ^ term range between 

magnitudes of about 1 and 5, whereas, the rise in the H20 band emittance 
occurs at magnitudes between 5 and 8 for reasonable values of mixing ratio 
in the stratosphere. 


In the design of an instrument system for mixing ratio inferences, 
these two curves should increase in approximately the same range of 


Log 


I 


term magnitudes. The implication of these results is plainly 
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seen when examining observed horizon radiances. The rise in the CO 2 pro- 
file occurs in the 30-km tangent height range, and the increase in the 
H 2 O radiances occurs at about 10 km (see figure 1) . Ideally the increase 
in the two radiance profiles should coincide in the same range of tangent 
heights. This may be accomplished by shifting the CO 2 band-pass filter 
to shorter wavenumbers (longer wavelengths). However, shifting the band- 
pass filter to more transparent regions of the CO 2 band entails some 
sacrifice of the signal-to-noise ratio at higher altitudes. 


Summary Remarks 

Analysis of the altitude dispersion characteristics of the horizon 
radiance in the CO 2 band indicates that the instant inversion technique 
will work for tangent heights greater than 20 km, and the accuracy should 
improve with altitude. It was shown that the emittance properties of the 
atmosphere appear to be linearly correlated with parameters characteristic 
of different air mass types. This property allows a calibration of the 
emittance for different climatological atmospheres and tangent height al- 
titudes. Deviations from the condition of thermodynamic equilibrium are 
quite small and are of the order of 1 percent for horizon radiances at 
60 km. Finally, it is shown from a comparison of CO 2 and H 2 O band emit- 
tances that the spectral intervals of scanner observations are not matched 
for purposes of moisture inferences using both bands. However, the fea- 
sibility of the concept is demonstrated by noting that the shapes of the 
emittance curves are similar. What is needed is a radiometer designed to 
sense CO 2 emission in more transparent regions of the band, i.e., at 
shorter wavenumbers or longer wavelengths. 
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RESULTS 


The instant inversion technique was applied to numerous horizon radi- 
ance profiles in the CO 2 band that were computed numerically on a digital 
computer- These test profiles were selected from NASA-CR66184 and NASA- 
CR6618'5 (ref. 13) so that the samples would represent the atmospheric vari- 
ations over the full range of seasonal and geographical conditions. The 
technique was used to infer temperature structure in the 20 to 60 km tan- 
gent height range, and the resulting temperatures were compared with ac- 
tual profiles used to compute the radiances. Statistics were compiled on 
the systematic error and random error uncertainties of several profile 
sets. These results allow a detailed evaluation of the instant inversion 
technique for inferring temperature structure. 

The technique was then applied to two radiance profiles in the CO 2 
band based on Scanner observations in order to illustrate its application 
to real data. 

A detailed analysis of H 2 O radiance profiles from Scanner observa- 
tions was performed to estimate average radiances as a function of lati- 
tude. These average radiances were then inverted to determine moisture 
profiles at different latitudes within the geographical area of observa- 
tion. An investigation was performed to determine the sensitivity of in- 
ferred water vapor mixing ratios to different transmittance models. In- 
ferred moisture profiles include the tangent height range from 10 to 35 
km. 


Technique Evaluation 

Criteria . - The instant inversion technique was applied to several 
independent sets of theoretical radiance profiles in the 615 to 715 cm“^- 
CO 2 band. Each selected set is noted by the following title: 

(1) 30 GCA Profiles - selected from reference 13 by the investi- 
gators . 

(2) 44 Aux. Profiles - the last 44 auxiliary profiles contained in 
NASA-CR66185 (ref. 13) which represent a collection of computed radiance 
profiles based on observed temperature and pressure profiles up to alti- 
tudes of 80 to 90 km. 

(3) 36 Aux. Profiles - a selected sample of 36 profiles contained 
in the 44 profiles of (2). 

(4) 30 LRC Profiles - a sample of 30 profiles provided by the scien- 
tists at the Langley Research Center. The origin of these profiles is un- 
known to the investigators of this report. 

Statistics were compiled on the accuracy of the inferred temperature 
profiles for each sample set. Both the systematic error se(h) and random 
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error re(h) were determined for nine tangent height levels from 20 to 60 
km. An entire set of calculations were performed for each of two emit- 
tance calibration parameters based on (1) the 10-mbar height and (2) the 
sum of nine radiances between the 20 and 60 km tangent height levels in 
increments of 5 km. Finally, statistical calculations were performed for 
three criteria of temperature comparisons; a comparison between the com- 
puted temperature and (1) the observed temperature at the 50-percent radi- 
ance level, (2) the average temperature over a 2-km layer centered about 
the level in question and (3) the average temperature over a 5 -km layer 
centered about the layer in question. 

Tables III and IV present tabulations of statistical results in terms 
of systematic and random temperature errors for 5-km layer comparisons. 

One obvious result from an inspection of systematic errors is a consis- 
tent underestimate of the temperature structure in the vicinity of the 
stratopause between 50 and 55 km. The error is on the order of 10° to 
13°C too low. 

The inability of the technique to predict temperature extremes fol- 
lows from the approximations made in the theory. The assumption is made 
that the variation of the Planck radiance above a tangent height may be 
approximated by an average Planck radiance B^ v [T(z')]* Thus, in regions 
of the atmosphere where there are temperature maxima, the technique will 
consistently underestimate the temperature. Conversely, in regions of 
temperature minima, the technique will overestimate the temperature. 
Similarly, but to a lesser extent, in regions of the atmosphere where 
the temperature is increasing above the tangent height level, the tech- 
nique will overestimate the temperature. Conversely, in atmospheric re- 
gions having a temperature lapse, the technique will underestimate the 
temperature. 

A practical solution to the problem of overestimating or underesti- 
mating the temperature is to simply correct for the systematic error. 

This adjustment would be based on statistical averages of the required 
corrections for many radiance profiles. 

It is noted that the systematic errors for 44 Aux. Profiles in Table 
III are all positive, considerably different than similar results for the 
30 GCA Profiles. Each of the 44 profiles in this set was examined in de- 
tail to ascertain the reason for the difference. It was determined that 
eight of the 44 profiles were associated with arctic air masses having 
extremely cold temperatures in the troposphere and warm temperatures near 
the stratopause. These profiles were significantly different than any of 
the 20 climatological atmospheres used to calibrate the emittance as a 
function of air mass. Thus the computed emittances were too low in mag- 
nitude. Had these profiles been included in the calibration analysis, 
the emittances would have been approximated by a quadratic polynomial in- 
stead of a straight line as in equation (15). 

After removing the eight profiles associated with arctic air masses, 
statistics were compiled for the remaining 36 Aux. Profiles. The resulting 
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TABLE III. - STATISTICAL EVALUATION OF INSTANT INVERSION TECHNIQUE 
Emittance Parameter - 10 mbar Surface 


Tangent 
Height h, km 

30 GCA 
se(h) 

Profiles 

re(h) 

44 Aux. 
se(h) 

Profiles 

re(h) 

36 Aux. 
se(h) 

Profiles 
re (h) 

30 LRC 
se(h) 

Profiles 
re (h) 

20 

* 

1.9 

* 

3.2 

2.4 

4.2 

3.6 

3.6 

1.5 

2.2 

25 

2.8 

2.3 

3.6 

3.3 

4.4 

2.7 

1.4 

1.9 

30 

2.1 

3.4 

2.3 

3.0 

2.7 

2.6 

1.0 

1.6 

35 

0.8 

1.5 

0.5 

3.8 

-0.7 

2.8 

-1.6 

2.3 

40 

-2.1 

2.8 

1.7 

12.6 

-3.7 

4.6 

-5.4 

6.4 

45 

-5.5 

4.9 

3.3 

25.8 

-7.6 

6.8 

-7.2 

4.7 

50 

-9.5 

6.5 

5.3 

38.7 

-10.6 

8.3 

-10.5 

5.4 

55 

-12.3 

7.2 

7.4 

52.1 

-12.8 

9.4 

-10.5 

5.6 

60 

-14.2 

6.2 

0.1 

35.8 

-14.4 

8.6 

-11.1 

5.8 


* 


Units of 



se(h) 


systematic error; re(h) 


random error. 


TABLE IV. - STATISTICAL EVALUATION OF INSTANT INVERSION TECHNIQUE 
Emittance Parameter - Summed Radiances 


Tangent 
Height h, km 

30 GCA 
se(h) 

Profiles 
re (h) 

44 Aux. 
se(h) 

Profiles 
re (h) 

36 Aux. 
se(h) 

Profiles 
re (h) 

30 LRC 
se(h) 

Profiles 
re (h) 

20 

* 

2.0 

* 

3.3 

2.2 

4.9 

3.7 

3.9 

2.0 

2.3 

25 

3.0 

2.6 

2.9 

5.2 

4.8 

3.4 

2.7 

2.2 

30 

2.5 

3.7 

0.6 

6.9 

3.2 

4.3 

3.4 

1.5 

35 

1.3 

3.1 

-3.0 

7.6 

0.1 

4.1 

2.7 

2.0 

40 

-1.4 

3.5 

-5.1 

7.0 

-2.5 

4.3 

1.1 

5.5 

45 

-5.1 

3.4 

-7.9 

5.7 

-6.0 

4.3 

1.4 

1.9 

50 

-8.4 

2.7 

-10.4 

5.1 

-8.8 

3.8 

-0.6 

2.1 

55 

-11.3 

2.3 

-12.5 

4.5 

-11.1 

3.4 

-0.2 

3.1 

60 

-13.3 

2.0 

-14.0 

3.7 

-13.0 

2.6 

-1.9 

4.9 


* 


Units of 



se(h) 


systematic error; re (h) 


random error. 


systematic errors are in excellent agreement with those for 30 GCA Pro- 
files. 

The systematic errors associated with the 30 LRC Profiles are differ- 
ent than the results of the other profile sets. This is especially true 
of the values presented in Table IV where the errors are slightly nega- 
tive in the 50 and 60 km tangent height range. The key parameter to de- 
termine in the inversion solution is the total emittance at each tangent 
height. The results indicate that the profiles in this group are signi- 
ficantly different than those considered in the emittance calibration 
statistics. The problem encountered with these profiles is symptomatic 
of the eight arctic air masses in the 44 Aux. Profile set. 

The inaccuracy of the emittance calibration is a major shortcoming 
of the instant inversion technique. One possible improvement that could 
be made is to develop a more sophisticated method of emittance calibra- 
tion. Probably a more fruitful avenue to pursue is the development of 
an inversion technique which converges to a solution by the method of 
successive approximations. In this case, the emittance would be recom- 
puted after each iteration until a solution is determined within the 
bounds of the instrumental error. An approach such as this would cer- 
tainly overcome the inherent shortcomings of the present technique. 

The random temperature error is a better indication of the prediction 
accuracy of the instant inversion technique, since it is possible to cor- 
rect for systematic errors. Results in Table III based on the 10-mbar 
surface indicate random errors of about 2 to 3°C for the 20 to 40 km al- 
titude range. This error is much larger in the vicinity of the strato- 
pause and may be attributed to a poor statistical coupling between the 
height of the 10-mbar surface and atmospheric structure above 50 km. 

Results in Table IV are more promising than those in Table III. Ran- 
dom errors have magnitudes on the order of 3 to 4°C for all altitudes 
tested, and there is an indication that the temperature accuracy improves 
for higher tangent heights. Even the results for the 44 Aux. Profiles are 
reasonable even though the aforementioned arctic air masses were included 
as part of the data sample. 

Statistical results were determined for three criteria of temperature 
comparisons - a comparison between computed temperatures and the observed 
temperature at the same level, and the average temperature over 2 and 5 km 
layers in the atmosphere. The change in random errors for the three cri- 
teria are presented in Table V. As can be seen from the results, the ran- 
dom error increases about 0.5°C when going from a 5-km layer to a point- 
by-point comparison. The fact that the instant inversion technique pre- 
dicts temperature as well as it does may be attributed to the geometry 
of limb-viewing and the peaked nature of the weighting functions. 
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TABLE V. - RANDOM ERRORS FOR THREE CRITERIA OF TEMPERATURE COMPARISONS 


Parameter - Summed Radiances 


Tangent 
Height h, km 

30 GCA Profiles 
Point 2 km 5 km 

36 Aux. Profiles 
Point 2 km 5 km 

30 LRC Profiles 
Point 2 km 5 km 


* 

V? 

* 







20 

4.4 

4.1 

3.3 

4.3 

4.2 

3.9 

2.7 

2.4 

2.2 

25 

2.8 

2.7 

2.6 

3.6 

3.5 

3.4 

2.4 

2.4 

2.1 

30 

4.4 

4.2 

3.7 

5.0 

4.8 

4.3 

1.7 

1.6 

1.4 

35 

3.5 

3.3 

3.1 

4.4 

4.3 

4.1 

2.1 

2.0 

1.8 

40 

3.8 

3.7 

3.5 

4.6 

4.5 

4.3 

2.2 

2.1 

1.6 

45 

4.2 

4.1 

3.4 

4.7 

4.6 

4.3 

1.9 

1.9 

1.5 

50 

2.7 

2.7 

2.7 

3.8 

3.8 

3.8 

1.3 

1.3 

1.5 

55 

2.4 

2.3 

2.3 

3.5 

3.4 

3.4 

2.4 

2.4 

2.4 

60 

2.0 

2.0 

2.0 

3.3 

3.0 

2.6 

4.3 

4.3 

4.2 


* o 

Units of C. 


Temperature Inferences from Scanner Observations 

The instant inversion technique was used to invert observations of 
horizon radiance in the CO 2 band in order to illustrate its application 
to realistic data. The data were extracted from a report (ref. 8) and 
are associated with average measured radiance profiles for Cells No. 4 
and 7. The radiances were inverted and then corrected for systematic 
errors according to the results in Table IV. The inferred temperatures 
are compared to those temperatures deduced from meteorological analyses 
for the Scanner Mission. The agreement is good, and the temperature 
error is in line with the magnitudes of random errors presented in Table 
IV. 


Moisture Inferences from Scanner Observations 

It was mentioned in the discussion of the results of figure 7 that 
the spectral intervals for CO 2 and H 2 O on the Scanner flight are not op- 
tically matched for purposes of inversion. Observations of the H 2 O hori- 
zon radiance profiles are available in the 10 to 35 km tangent height 
range. Unfortunately, temperature inferences using the instant inversion 
technique from the CO 2 band are not applicable below the 30-km level. 
Therefore, it is not possible to perform every detail indicated by the 
development of the inversion theory. 

In order to perform the inversion of H 2 O radiance profiles, it is 
necessary to know the emitting temperature and the ^ term in equation 

(13). Normally these parameters would follow as results from the inver- 
sion of CO 2 radiances. The magnitude of these terms in the mixing ratio 
computation [equation (14) ] were determined from available temperature 
and pressure data (ref. 8) associated with the observed profiles. The 
effective temperature was computed as an average, weighted by the den- 
sity along the line-of -sight , and the ^ term was computed using the 
formulation of equation (11). 

Analysis of Scanner Data . - The geographical coverage of the Scanner 
measurements is indicated by the hatched area of figure 10 - also shown 
are the nephanalysis and significant weather at time of observation. In 
reference 8 the radiance profiles were partitioned into seven geographi- 
cal cell locations which vary as a function of latitude along the geo- 
graphical annulus of observation. 

An analysis was performed during the course of this study to deter- 
mine average horizon radiances for the 315 to 475 cm"^ H 2 O band as a 
function of latitude. Average radiance observations for this band were 
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ALTTTUDE z,km 



TEMPERATURE T, °K 


Figure 8. Comparison of inferred temperatures with the observed profile 
for Cell No. 4 of Scanner data. 


30 



ALTITUDE z, km 



Figure 9. Comparison of inferred temperatures with the observed profile 
for Cell No. 7 of Scanner data. 
















plotted in figure 11 for six tangent height levels from 10 to 35 km. 

Trend lines drawn through the scatter of points, and radiance magnitudes 
were determined for the latitude nosition associated with each cell. Re- 
sulting average measured radiance profiles are tabulated in Table VI. 

The average radiances for the 10 km tangent height should be viewed 
with some caution. The discussion in reference 8 indicates that both 
high level clouds and instrument response could introduce systematic biases 
into the magnitudes of the 10 km radiances. These bias factors have a tend- 
ency to lower the magnitude of observed radiances which in turn would result 
in decreased magnitudes of computed mixing ratios. Even with this qualifi- 
cation, however, the computed mixing ratios appear to have reasonable mag- 
nitudes . 

Effective emitting temperatures were derived from model atmospheres 
developed (ref. 8). These temperatures were then converted to Planck radi- 
ances [equation (2) ] representative of the 160 cm"^ band-pass at a central 
wavenumber of about 370 cm‘^. This wavenumber is shifted 25 cm“^ off the 
center of the band-pass filter to account for the asymmetry of the H 2 O 
band within the 315 to 475 cnT^ spectral interval. Results are presented 
in Table VII. 

Comparison of transmittance models . - An investigation was performed 
to determine the sensitivity of the inferred moisture profiles to differ- 
ent transmittance models. This analysis was based on a comparison of 
Elsasser's model (ref. 12) and that of Smith (ref. 14). Calculations of 
H 2 O band emittance were made for a variety of mixing ratio profiles using 
the Standard Atmosphere temperature profile for both transmission models. 

The results are presented in the graphs of figure 12 where the band emit- 
tance is plotted as a function of [equation (12)] on log-log scales. 

It is noted that Smith's model yields smaller values for emittances 
greater than 0.06 and larger values for emittances smaller than 0.06 in 
magnitude. An emittance of 0.06 corresponds to a tangent height of roughly 
25 km. The emittance model can have a pronounced effect on the inference 
of moisture profiles in both the lower and upper stratosphere as will be 
shown in the next section. Smith's emittance model was used to compute 
moisture profiles for this study since it is based on recent calculations 
and has a finer spectral resolution than Elsasser's model. 

Inferred moisture profiles . - Calculations of the water vapor mixing 
ratio were performed using the data in Tables VI and VII, Smith's trans- 
mittance model in figure 12 and equation (14) . The results are presented 
in Table VIII for the seven cells. An average of the mixing ratios is 
shown in the last column. Computed mixing ratios indicate a general de- 
crease to about 0.006 g/kg in magnitude at the 23-km level, a gradual in- 
crease to about 0.009 g/kg at 32 km and then a decided decrease to about 
0.005 g/kg at 37 km. 

The results in Table VIII are plotted in terms of latitudinal vari- 
ations in figure 13. Some cautious conclusions may be drawn from these 
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TABLE VI. “ AVERAGE MEASURED RADIANCE PROFILES 
FOR 315 CM" 1 TO 475 CM" 1 


Tangent 
Height, km 


Radiance N H (h) 

, watts/m^-sr (eq. 

• 8) 


Cell 1 
58°N 

Cell 2 
53°N 

Cell 3 
47 °N 

Cell 4 
43 °N 

Cell 5 
35°N 

Cell 6 
21°N 

Cell 7 
17 °N 

10 

6.20 

6.50 

6.70 

6.60 

6.90 

7.10 

7.30 

15 

2.80 

2.80 

2.90 

3.10 

3.30 

2.90 

2.80 

20 

1.74 

1.64 

1.59 

1.63 

1.60 

1.42 

1.40 

25 

0.90 

0.90 

0.89 

0.97 

0.93 

0.79 

0.80 

30 

0.51 

0.46 

0.55 

0.59 

0.55 

0.46 

0.48 

35 

0.28 

0.23 

0.25 

0.27 

0.29 

0.26 

0.25 


35 




TABLE VII. - EFFECTIVE EMITTING TEMPERATURES AND PLANCK RADIANCES 
FOR SCANNER OBSERVATIONS FROM 315 CM" 1 TO 475 CM" 1 


Tangent 

Height, 

km 



Temperature, 

°K 


• 8) 

Planck 

Radiance 

B h [T(z')] 

, watts/m 4 

2 (eq< 

-sr 

Cell 1 
58°N 

Cell 2 
53°N 

Cell 3 
47 °N 

Cell 4 
43°N 

Cell 5 
35°N 

Cell 6 
21°N 

Cell 7 
17°N 

i n 

224.0 

223.1 

223.4 

219.6 

220.9 

223.1 

223.1 

1U 

9.89 

9.78 

9.82 

9.27 

9.53 

9.78 

9.78 


225.5 

224.2 

220.8 

218.6 

215.7 

207.8 

207.7 

Id 

10.07 

9.92 

9.52 

9.26 

8.96 

8.08 

8.07 

on 

228.2 

227.4 

226.6 

224.5 

222.4 

219.1 

218.9 


10.39 

10.30 

10.20 

9.95 

9.70 

9.32 

9.30 


234.6 

231.2 

231.8 

230.9 

230.2 

229.4 

229.0 

Z.D 

11.14 

10.74 

10.81 

10.70 

10.63 

10.53 

10.48 

o n 

240.9 

240.3 

241.1 

240.9 

239.1 

237.1 

235.6 

JU 

11.93 

11.86 

11.95 

11.93 

11.70 

11.46 

11.26 


250.6 

249.3 

249.6 

249.9 

251.5 

250.6 

250.4 

35 

13.14 

13 . 00 

13.04 

13.08 

13.24 

13.14 

13.12 


36 









TABLE VIII. - INFERRED WATER VAPOR MIXING RATIOS FROM SCANNER OBSERVATIONS 


Effective 

Altitude, 

km 



Mixing 

ratio Wjj, 

g/kg 



Mean 

Values 

Cell 1 
58°N 

Cell 2 
53°N 

Cell 3 
47 °N 

Cell 4 
43°N 

Cell 5 
35°N 

Cell 6 
21°N 

Cell 7 
17°N 

14.9 

0.0105 

0.0145 

0.0166 

0.0178 

0.0204 

0.0195 

0.0224 

0.0174 

18.6 

0.0044 

0.0046 

0.0055 

0.0071 

0.0095 

0.0089 

0.0081 

0.0069 

22.8 

0.006 

0.0059 

0.0058 

0.0068 

0.0069 

0.0060 

0.0059 

0.0062 

27.2 

0.0068 

0.0076 

0.0071 

0.0087 

0.0066 

0.0066 

0.0068 

0.0072 

31.8 

0.0076 

0.0063 

0.0089 

0.0112 

0.0105 

0.0078 

0.0089 

0.0087 

36.6 

0.0051 

0.0034 

0.0036 

0.0054 

0.0063 

0.0050 

0.0042 

0.0047 









trend curves. Latitudes north of 45 degrees are consistently drier than 
those to the south. The subtropical region is quite moist up to about 
20 km; however, the zone between 35 and 45°N is wettest above the 20-km 
level. It is interesting to note that this region is associated with a 
frontal weather situation in the troposphere (see figure 10). One may 
speculate that the weather activity has injected moisture into the 
stratosphere . 

Figure 14 shows a comparison of mean mixing ratio profiles for two 
transmittance models. The results of Table VIII were based on Smiths 
transmittance model - similar calculations were performed in order to de- 
termine a mean profile for Elsasser's model. It is interesting to compare 
these profiles with the emittance curves in figure 12. The two profiles 
cross at an altitude of 27 km, which corresponds to the same intersection 
of the emittance curves. Above this level, larger emittance values yield 
small magnitudes of mixing ratio, and the converse is true for altitudes 
below 27 km. 

The mixing ratios derived by this technique are generally larger in 
magnitude than typical values of between 0.002 and 0.003 g/kg observed by 
Mastenbrook (ref. 15) for the 25 to 30 km altitude range. However, some 
of his observations are considerably larger than the magnitudes derived 
here. It is difficult to establish at this time whether the inferred 
moisture profiles are too large in magnitude since independent measure- 
ments are not available at the time of observation. 


Summary Remarks 

The instant inversion technique was evaluated statistically using 
three independent sets of radiance profiles. Results indicated that the 
technique consistently underestimated temperature maxima and overesti- 
mated temperature minima. However, it is possible to correct inferred 
temperatures for systematic errors by adding or subtracting a given mag- 
nitude of temperature based on an average for all sets of profiles. Ran- 
dom temperature errors were of the order of 3 to 4°C for each set of pro- 
files. The random error will increase to about 5°C if the uncertainty in 
the systematic error correction is Included in the calculations. The in- 
stant inversion technique works quite well in view of its simplicity of 
approach. The technique was successfully applied to observed horizon 
radiances from Project Scanner in both the CO 2 and H 2 O rotational bands. 
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ALTITUDE 



Figure 14. Comparison of mean mixing ratio profiles inferred from horizon 
radiances measurements 315 to 475 cm - ^ (ref. 8) for two trans- 
mission models. 
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SUMMARY ANALYSES 


Part of the effort under this program was directed toward an 
analysis of errors pertinent to the inference of temperature and mois- 
ture from horizon radiance measurements. This task was approached by 
considering representative errors of radiometer systems that could be 
flown on earth-orbiting spacecraft. Calculations were based on the 
formulation specified by the instant inversion technique. Random error 
uncertainties were estimated for both temperature and moisture infer- 
ences as a function of altitude in the atmosphere. 

Analyses were also performed to determine optimum spectral inter- 
vals in the and CO^ rotational bands for radiometer systems to be 

used for purposes of temperature and moisture inferences. 

Results of these analyses are presented in this section of the 
report. 


Achievable Temperature Inference Accuracies 

It has been demonstrated that the instant inversion technique is 
a reasonably good predictor of temperature and moisture profiles. Since 
the formulation is simple, it is quite useful in performing error analy- 
ses of prediction accuracies. There are two sources of observational 
errors — radiance measurement errors and uncertainties in specifying 
the tangent height altitude. There are also errors associated with the 
emittance e(h) or transmission model used in the calculations and errors 
in the inversion technique itself. It is reasonable to neglect errors 
due to the inversion scheme since an iteration approach with regression 
analyses will certainly overcome the shortcomings of the instant inver- 
sion technique. Differentiating equation (5) gives 


Assuming these errors are uncorrelated, the rms temperature error is 
expressed by 


AT(h) 


rms 



1 

e(h) 


AN (h) 2 





1/2 

( 18 ) 


A random radiance accuracy of 0.01 watts /m^ sr and a tangent height error 
of = 0.4 km are assumed to be reasonable for instrumentation on space- 
craft. The random uncertainty in the emittance e(h) was estimated to be 


about 2 percent 



0.02 




These values along with total emittances 
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and radiances . for the Standard Atmosphere in the 615 cm“l to 715 cm"^ 
spectral range were substituted in equation (18) to obtain estimates of 
the rms temperature error. The results are presented in figure 15. 

Errors are generally 1°C in the lower stratosphere, increase gradu- 
ally to about 4°C near the stratopause and further to about 5°C at the 
60-km level. The accuracy of the tangent height level is the principal 
contributor to the rms temperature error between 20 and 55 km, and the 
measured radiance accuracy becomes the dominant source of error above 
the 60-km level. These estimated rms temperature accuracies are cer- 
tainly competitive with those from meteorological rocket soundings. 
Wagner (ref. 16) indicates temperature errors on the order of 5°C to 
9°C in the 50 to 60-km altitude range for bead thermistor sensors used 
by the Meteorological Rocket Network. 

It may be argued that the determination of the tangent height is 
more systematic than random in character for an individual profile. 

This follows from an argument that the spacecraft's attitude may be un- 
certain by about 0.4 km; however, this uncertainty may be constant over 
the time interval of measurement (less than 1 sec). If this situation 
is realistic, the entire radiance profile is in error by a systematic 
increment of tangent height. In a corresponding manner, the inferred 
temperature profile may be shifted by a similar error in the altitude 
from of reference. On this basis the temperature gradients with alti- 
tude may be quite accurate, but the absolute temperature at a given 
level will be in error by magnitudes indicated in the analysis. 

Possible verification of this systematic error and its effect on 
temperature accuracies could be accomplished in a more detailed analysis 
of the problem. 


Achievable Moisture Inference Accuracies 


Estimates of the prediction accuracies of water vapor mixing ratio 
may be determined by an approach similar to that used for temperature 
accuracies. From equations (8) and (10), 


N„(h) 

t (h) = 1 - ^ (19) 

H B (h) 


Differentiating equation (19) gives 


dT R (h) 


d V h > , 1 

B R (h) B H (h) 


dN R (h) 

Sh 


N (h) 

dh + zr o dB (h) 

V h) 


( 20 ) 
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Assuming the transmittance is approximated by the following exponential 
form 


V h > - 



( 21 ) 


differentiation gives 


dT R (h) 



T R (h) dW R 


( 22 ) 


From equation (21) 



(23) 


Combining equations (22) and (23) 


dW R dT R (h) 

W “ " V h > (^w) 


(24) 


Assuming the errors are uncorrelated, the fractional rms mixing ratio 
error may be expressed from equations (20) and (24) by 



Equation (25) was evaluated over the altitude range from 10 to 50 km 
using realistic magnitudes of errors and input values. A random radiance 
accuracy of 0.025 watts/m^ sr and a tangent height error of 0*4 km were 
assumed. Values for the transmittance Tjj(h), horizon radiances Nu(h) and 
Planck radiances B^(h) were extracted from the scanner data in tables VI 
and VII. Theoretical calculations of these parameters were used for data 
in the 35 to 50-km altitude range. Temperature errors (or AB^(h) errors) 
were based on data in figure 15. The results of the calculations are pre- 
sented in figure 16 as percent errors of the absolute values of the mixing 
ratios. 
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RMS mixing ratio errors are generally 5 to 8 percent in the 10 to 
30 km altitude range. The percent error increases consistently above the 
30-km level to a value of 30 percent at 51 km. In order to obtain the 
absolute error from these result s, the value of the mixing ratio is mul- 
tiplied by the percentage for a given level. For example, a 20 percent 
error for a mixing ratio of 0.010 g/kg amounts to a rms uncertainty of 
0.002 g/kg. If the mixing ratio were 0.002 g/kg, the rms error would be 
only 0.0004 g/kg. 

The calculations indicated that the tangent height uncertainty in 
the 10 to 25 km altitude range was the prime source of error. This is 
the region of the ^0 radiance profile where the slope of the curve 
(dN H /dh) is greatest in absolute magnitude (see figure 1). Above 25 km, 
the radiance accuracy is the dominant source of error. The calculations 
also indicated that errors due to temperatures are relatively unimportant 
compared to radiance and tangent height errors. This means that the mix- 
ing ratio calculation is insensitive to temperature errors determined from 
CC >2 radiance profiles, and this result emphasizes the feasibility of using 
a dual radiometer system to obtain temperature and moisture profiles in 
the stratosphere. 


Optimum Spectral Intervals 


The general criterion that may be applied to the selection of band- 
pass filters is the maximum radiant power that could be intercepted by 
a radiometer (i.e., maximum signal -to -noise) . A complete treatment of 
this problem would consider a detailed analysis of entire radiometer sys- 
tems. The approach to the problem here is restricted in scope and con- 
siders only atmospheric physics, and thus disregards instrumental prob- 
lems such as filter response and detector response. 


CO? band . — Two factors influence the problem of selecting CO 2 band- 
pass filters: (1) maximum possible signal for measurements at high tan- 

gent height levels, and (2) an optical matching of the radiance profile 
associated with the H^O band (see section on a Comparison of CO 2 and ^0 
Band Emi ttances) . The investigation approached the problem by consider- 
ing optical path lengths as a function of tangent height, and then re- 
lating these to the optical properties of the band. It is also desirable 
to maintain the peaked nature of the weighting functions at all altitudes 
(see figure 3), which maximizes the inversion possibilities. These con- 
siderations led to the following suggestions for the selection of band- 
pass filters 


(1) Single filter radiometer — spectral response 590 cm‘^ to 700 cm"^, 

(2) Dual filter radiometer — spectral responses 615 cm"^ to 700 cm“l 
and 570 cm" to 660 cm“^ . 


The selection of the single filter is a compromise and attempts to fulfill 
the two requirements of maximum signal and optical matching. 
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The selection of two filters is quite attractive in that each filter 
optimizes a given altitude region of the atmosphere. The 615 cm“l to 700 
cm~l filter applies to the 30 to 65-km altitude range, and the 570 cm~^ 
to 660 curl filter provides good signals in the 10 to 40-km altitude range 
The overlap in the 30 to 40-km range allows a continuity of measurements 
for purposes of temperature inferences . 

It should be noted that the spectral region beyond 700 cm"^ is 
avoided due to the presence of a moderate ozone absorption band at 14 mi- 
crons. Secondly, there is justification for extending the spectral re- 
sponse towards shorter wavenumbers in the wings of the H 2 O band since 
there are relatively small path lengths of water vapor in the stratosphere 
even when considering limb viewing. Reduced path lengths at 10 km are gen 
erally about 0.1 pr cm, an order of magnitude less than the total preci- 
pitable water concentration of a vertical column in the troposphere. Cal- 
culations of the radiance from water vapor in the 570 to 600 cm”! spectral 
range indicate that the contribution by H 2 O to the total radiance is on 
the order of 1 percent at 10-km tangent heights. Therefore, it is reason- 
able to extend the band-pass filter to 570 cm~l in order to obtain maximum 
signal at altitudes of 10 to 15 km. 

H? 0 band . - Several factors can influence the selection of a band- 
pass filter in the H 2 O rotational band. The principal factor is the op- 
tical properties of the band as a function of spectral region. Other 
factors such as effective emitting temperature, tangent height and width 
of the band-pass filter also influence the selection. It is emphasized 
again that the investigation here considers only the atmospheric physics 
part of the problem and disregards the instrument itself. 

The investigation was approached by considering the detailed spectral 
emitting properties of the band from 200 to 600 cm"^ for a number of dif- 
ferent atmospheres and water vapor concentrations. Calculations were 
based on equation (8), using Smith’s transmittance model (ref. 14) to com- 
pute the emittance. The computed spectral radiances were smoothed by dif- 
ferent band-pass filter widths of Av = 40, 80, 120, and 160 cm“^ in a man- 
ner which simulates what the radiometer would receive. These results were 
then analyzed to determine the optimum filter configuration as a function 
of tangent height. The study is limited to filters positioned in the 20p 
to 3 5p. region, in accordance with the statement of work. 

The calculations indicated that the spectral radiance over spectral 
widths Av of 5 cm~l exhibits considerable fine scale structure as is in- 
dicated by the results in figure 17 for a 20-km tangent height. Spectral 
radiances like those in figure 17 were computed for tangent heights from 
the surface to 50 km in 5 km intervals, and for a variety of selected at- 
mospheric conditions. These detailed radiances were smoothed, and figure 
18 is an example of results when the radiances were smoothed with the 
four filter functions. The selection of the optimum central wavenumber 
for a band-pass filter is based on the region which gives the maximum 
radiant power. 
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SPECTRAL RADIANCE, Watts / m‘ sr 





The final results of this analysis are presented in figure 19 
which indicates the filter position as a function of tangent height for 
the H 2 O rotational band. For all practical purposes, the filter posi- 
tion is independent of band -pass width. The controlling factor in the 
altitude range from the surface to 15 km is the effective emitting 
temperature which manifests itself in terms of the nonlinear Planck 
radiance. This effect applies only in regions of the atmosphere where 
the optical path length is nearly opaque. As opacity decreases from 
near unity, the position of the filter rapidly moves to longer wave- 
lengths. This characteristic is a result of the optical properties 
of the band . 

It is concluded from this analysis that the band-pass filter should 
be positioned at 35\i or even longer wavelengths, if possible. This 
filter position provides the maximum radiant power to the radiometer 
and enhances the accuracy of inferred moisture profiles. 
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TANGENT HEIGHT 






CONCLUDING REMARKS 


A detailed investigation has been performed to assess the feasi- 
bility of inferring stratospheric temperature and moisture profiles from 
horizon radiance measurements. Both the 15ji CO 2 band and the H 2 O rota- 
tional band were considered as potential spectral regions in the infrared 
whose radiance measurements can be used for temperature and moisture 
inferences. It was concluded that useful temperature measurements can 
be obtained from CO 2 band radiances over the altitude range from 10 to 
65 km. Water vapor mixing ratios can be determined from H 2 O band radi- 
ances for altitudes from about 10 to 50 km. 


The optical properties of the CO 2 band (615 to 715 cm"^) and the 
H 2 O band (315 to 475 crn'l) were studied in detail. An inference method 
termed the instant inversion technique was developed for both bands and 
then applied to theoretical profiles as well as observed radiance pro- 
files from Project Scanner. There was good agreement between computed 
and observed temperature profiles, and the inferred magnitudes of water 
vapor mixing ratio were in general agreement with recent observations 
from balloons by frost-point hygrometers. 


Statistical analyses of temperature inference results for approx- 
imately 100 sample profiles indicate that the instant inversion technique 
can predict temperatures to accuracies of about 5°C. The error is due 
principally to the estimation of total atmospheric emittance at a tangent 
height. This shortcoming may be eliminated by developing an inversion 
technique which converges to a solution by the method of successive 
approximations. The emittance would be recomputed after each iteration 
until a solution is determined within the bounds of the instrumental 
error. A similar iterative approach may be used to infer moisture pro- 
files once the temperature profile is determined. 


Analyses were performed to estimate temperature inference accuracies 
for realistic radiance profile measurements. Sources of random error 
included radiance measurements, tangent height error and uncertainties 
in the transmission model. Results indicate rms accuracies of about 1°C 
in the lower stratosphere, 4°C near the stratopause and 5°C at the 60-km 
level. The principal source of error was due to the tangent height 
determination. It may be argued that this error is more systematic than 
random in character for a given radiance profile which tends to reduce 
the random temperature error. Possible verification of this systematic 
error and its effect on temperature accuracies could be accomplished in 
a more detailed analysis of the problem. 

A similar study was carried out to determine moisture inference 
accuracies. Sources of error included radiance measurements, tangent 
height error and uncertainties in the temperature inference. Results 


53 



indicate rms accuracies on the order of 5 to 8 percent in the 10 to 30 km 
altitude range, and the error increases to about 30 percent at the 50 km 
level. The absolute error is interpreted to be the product of the per- 
centage and the absolute value of the mixing ratio. The tangent height 
error contributes to the inaccuracies in the lower altitude range, and 
the radiance error is dominant above 30 km. The calculations indicated 
that moisture estimates were insensitive to temperature errors. 

Analyses were performed to determine the optimum spectral intervals 
in the infrared whose radiance measurements can be used to infer temper- 
ature and moisture profiles in the stratosphere. It was suggested that 
two band -pass filters be used when making the CO 2 radiance measurements, 
one with a spectral response 615 to 700 cm~l and the second with a re- 
sponse of 570 to 660 cm’l. These filters would apply to the upper and 
lower stratosphere, respectively. A single broad band -pass filter is 
required for the H 2 O rotational band, and it should be centered at 35|j 
or even longer wavelengths, if possible. 

Horizon radiance observations offer an attractive means to sound the 
stratosphere. Global measurements of stratospheric temperature and 
moisture profiles are not currently available for forecasting or research 
use. The results of this study have demonstrated the feasibility of 
making these measurements from earth-orbiting spacecraft. Furthermore, 
it should be possible to determine profiles of atmospheric ozone based 
on horizon radiance measurements in the 9.6j a band. Further studies 
should be directed toward the investigation of optimum inference tech- 
niques and specification of measurement requirements of instruments 
with the goal of developing the experiment for application on earth- 
orbiting spacecraft. 
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APPENDIX A 


TEMPERATURE INFERENCE AND EVALUATION COMPUTER PROGRAM 


A computer program termed INVERT was developed during the course 
of the study effort and was used extensively to test and evaluate the 
instant inversion technique. The following four points described the 
program’s capabilities and options. 

(1) Evaluates up to 100 radiance profiles at once. 

(2) Performs the instant inversion technique for 

(a) two spectral intervals 615-715 cm"^ and 600-725 cnT^, 

(b) two emittance calibration parameters, 

(c) two types of effective emitting levels. 

(3) Computes temperature errors for the above combinations of 
parameters and for a 

(a) point temperature comparison, 

(b) 2 -km layer average temperature comparison, 

(c) 5-km layer average temperature comparison. 

(4) Compiles statistics of temperature errors for all combinations 
of parameters and for the following statistics: 

(a) the root-mean-square (rms) error, 

(b) the systematic error component of the rms error, 

(c) the random error component of the rms error. 

INVERT is written in conventional FORTRAN II language and should 
operate without modification on most digital computers. The program 
is available for anyone desiring its use. Requests for source deck 
listing and operating instructions should be directed to the authors 
or to the technical monitor at LRC, Mr. R. E. Davis. 
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